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Abstract 

Using high-quality gas phase electron scattering calculations and multiple scattering theory, 
we attempt to gain insights on the radiation damage to DNA induced by secondary low-energy 
electrons in the condensed phase, and to bridge the existing gap with the gas phase theory and 
experiments. The origin of different resonant features (arising from single molecules or diffraction) 
is discussed and the calculations are compared to existing experiments in thin films. 

PACS numbers: 
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I. INTRODUCTION 



Since it was discovered^ that low energy electrons (LEE) can cause strand breaks in 
the DNA duplex, the interest in electron interactions with DNA, which originates preva- 
lently from its importance in radiation damage to living tissue and radiotherapy, has grown 
consistently-. The impact of ionizing electromagnetic radiation on matter causes the emis- 
sion of highly energetic electrons; these latter ionize the medium via an electromagnetic type 
of interaction, thus producing secondary electrons in large numbers^. Most secondary elec- 
trons are created with low energy (E<20 eV) and have a distribution with a most probable 
energy of 9 eV— . If the electron energy is above the DNA ionization threshold (7-10 eV) these 
electrons can ionize it, while over the entire 0-15 eV range they can be captured in a resonant 
anionic state. Fast, efficient dissociation pathways are known to exist for organic molecules 5 
from electronic dissociative states, dissociative electron attachment, dipolar dissociation and 
dissociative ionization, which can be related to capture or transfer of an electron or hole 
into a dissociative state.- 1 ^ 

In the context of sub-ionization threshold electrons, the importance of resonances has be- 
come evident. Specifically, electron capture by a DNA subunit to form a molecular resonance 
enhances greatly the rupture of chemical bonds within the molecule either by dissociative 
electron attachment (DEA) or the decay of the transient anion into a dissociative electron- 
ically excited state. These phenomena are reflected in the measured yield of strand breaks 
which exhibit prominent resonant features as a function of energyS. 

In a condensed environment, LEE are created within or outside DNA. The former nec- 
essarily have a high probability to interact with DNA. LEE created outside DNA may also 
interact with this molecule depending on the timescale of thermalization of the electrons. 
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In this case, other types of damage, such as radical-mediated damage, can become impor- 
tant or even predominant. In general, if the thermalization is slow, DNA is impacted by a 
relatively hot electron (meaning in this context an electron not fully solvated) from outside; 
this entails usually a timescale between 0.1 and 100 fs.- Others created outside the DNA do 
not arrive at the target with any appreciable energy and therefore cannot directly damage 
it, but they can form radicals along the way and these in turn can attack the DNA—. The 
prevalence of one or the other of these mechanisms probably is determined by the DNA 
concentration: if it is high, then the electrons are more likely to be created within DNA 
or impact it before being completely solvated. In principle, modeling the dynamics of the 
radiation tracks^ can yield very useful informations in this regard. 

Since radiation damage in a cell is a complicated problem, researchers have tried to sim- 
plify its description by observing first what happens to DNA components in the gas phase 
when impacted by LEEs; this literature is by now rather vast. The experimental approaches 
have been focused mostly on predicting products of dissociative electron attachment from 
DNA subunits and determining the role of shape resonances^^^i. The theoretical com- 
munity has provided information on resonances in elastic and inelastic scattering^^^^i22. 
While dissociative electron attachment calculations for these large targets are still beyond 
reach, quantum chemical methods have been employed to predict the weakest bonds in elec- 
tron attachmentri>2i and study the hydration effects on the anionic compounds 22 . It must be 
noted, however, how none of these methods will be able for the foreseeable future to deal 
with a molecule of tens of thousands of atoms like bacterial DNA in a process in which both 
electronic and nuclear degrees of freedom are involved, the scattering electron is unbound 
and the solvent has to be considered as well. While DFT is able to deal with large systems, 
nuclear dynamics on DFT surfaces is not as well defined 2 ^ 2 ^ and a scattering theory using 
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DFT is just starting to be developed 2 ^. 

This work aims to be a first step towards systematizing what happens in passing from the 
gas to the condensed phase. In particular, we will be concerned with the solid state, since 
this gives us the opportunity to relate to radiation damage experiments in thin film-' 2 and, 
equally importantly, to neglect fluctuations in the DNA structure, which are of fundamental 
importance in solution^ 2 ^. 

To this extent, the first step is to understand what happens inside the DNA polymer 
itself, which is represented here as rigid and immersed in vacuum. In this context it is 
possible to explore parameters like coherence length (which controls the exponential decay 
of the amplitude of the phase coherent part of the electron wave function), DNA sequence, 
and conformation, which could be important in radiobiological damage. 

We pursue this goal using a model, described below, which unites a recently developed 
multiple scattering framework— i2iL2ii2i with accurate electron scattering calculations for the 
DNA bases performed using the R-matrix methocV^ 2 -. 

In our model we have left out of the picture, for the moment, all the properties of the 
liquid phase, from fluctuations to the motion of the solvated electron in water. We plan to 
introduce the structural water (also present in thin films) as a scattering dopant in further 
studies using this model, which will be fairly straightforward. To introduce the liquid phase 
properties 33 would be instead much more difficult and require some essential modification 
to this framework. 

II. MODEL 

The multiple scattering framework we shall be using to get information on the elastic 
scattering of low-energy electrons (LEE) has been described in the series of articles by Caron 
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and Sanche^'SiLSi'^i. It was developed at a time when no scattering matrix information was 
available for the DNA bases. The objective of the toy model used was to find trends in cross 
sections and capture amplitudes resulting from regularity or disorder in the helix and base 
sequencing. With the advent of recently available scattering information on DNA bases, 
it is important to revisit elastic LEE scattering on DNA to get more precise numerical 
estimates and ascertain whether or not the previous conclusi ons are still valid. One obvious 
difference of the current calculations with respect to Refs. 



28 



29 



30 



31 



is the presence of 



shape resonances that carry over from the gas phase scattering of the subunits, which makes 
it possible to observe the fate of local resonances in the conjugate. We shall first study, in 
this paper, the idealized B-form of a GCGAATTGGC decamer (without backbone)^ and its 
regularly sequenced cousin, the poly (A) - poly (T) base pairs decamer. We shall then examine 
the A-form of the GCGAATTGGC decamer 34 . These structures are chosen to be ideal ones, 
although they are in general known to be slightly different^. Our prototype systems are 
decamers because 10 is the number of base pairs in one turn of the helix of B-form DNA. 
But let us first review the theoretical framework. 



A. Multiple Scattering Theory 



In Refs |28| and |29_|, we presented the basic equations for multiple electron scattering within 
macromolecules, including DNA. For the latter, we proposed a simple model of molecular 
subunits (i.e. bases, sugars, and phosphates) immersed in an optical potential U op , which is 
constant between their R- matrix shells (or between the muffin-tins), a working hypothesis 
that has been used in the calculations for simple molecules^, in the theory of low-energy 
electron diffraction (LEED) in solids^ and nanoscale structures^ . The only function of the 
real part of the optical potential is to account for the average energy seen by an electron. One 
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can quite generally describe the scattering problem of a molecular subunit by its scattering 
matrix Sl'l—^- where L = (/, m) are the angular momentum quantum numbers. Each 
molecular subunit has an incident plane wave of momentum k impinging on it plus the 
scattered waves of all other subunits. More specifically, we described the asymptotic form 
of the total wave function ^l n '(f) for a molecule centered at R n outside the R-matrix shell 
by the following equation 



^l n) (r) = 4rre 



LL' 



ji (kr n ) 5 UL + i (s^l - 6 l >l) hf 3 (kr r , 



where Yl are spherical harmonics, ji and hy are the spherical Bessel function and Hankel 
function of the first kind respectively, f n = f — R n , and 



4? = n (%) + \ E E 



kL 2 



2^2 



X (-1)' 



mi ,m,—m,2 ' 



(2) 



where 



mi ,m,— m' 2 



[4tt(2/ 1 + 1)(2/ + 1)(2^ + 1)] 







mi m — m^/ 



i? n — i? n /. Equation (T2J) implies a 



and (4m-4j is the Wi 

gner 3-j symbol^, and R nn 
coupled set of linear equations for all , which measure the resultant of the superposition 

fc L 

of the incident plane wave and the contribution from all other scatterers. As mentioned 
before^ 1 ^, the loss of coherence of the electrons due to inelastic collisions and perhaps also 
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to the presence of parasite scatterers (e.g. the water molecules in the grooves could be 
considered as such) can be invoked through an imaginary part in the background optical 
potential U 0] r^, i.e. an imaginary part to the electron wave number Im(A;) = Here £ 
acts coherence length for the electrons. 

B. Electron capture and scattering 

In an effort to extract physically meaningful information from the multiple scattering 
formalism, we had previously targeted a calculation of the capture amplitude vS n ' of an 
electron in a shape or core excited resonance of a basic sub unit positioned at R n . We had 
assumed a dominant capture channel symmetry corresponding to L and had used the one- 
center approximation of O'Malley and Taylor 42 . When generalized to a multiple scattering 
situation, this leads to 



where Vl is an energy and nuclear coordinate dependent amplitude. There is unfortunately 
no available theoretical information on the nuclear part of the wavefunction for the DNA 
bases at this time. So we shall only focus on the electronic part. 
We proposed 30 a weighted partial capture factor 




(3) 



(4) 



where the constituent partial capture factor 
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measures the partial wave decomposition of the total wave function averaged over the differ- 
ent bases positioned at R n . This would serve as meaningful measure of the effect of multiple 
scattering on the capture probability in the l a channel since j(l , R n ) = 1 for a lone plane 
wave. Note that 7(0, R n ) equals the average absolute square of the wave function at R n and 
thus 1^(0) measures the absolute square of the wave function averaged over all bases. 

The total elastic cross section, for a finite size macromolecule, would also be of interest. 
This is, however, a somewhat elusive quantity when there are losses. Technically, we can 
expand the scattered part of Eq. flTJ) around the geometric center R GC of the macromolecule. 
In this reference system, remembering that f n = r — R n , one has 



Yjj (fi,J hf (kr n ) = £ ^-^-ir'^X-rn^ (" 



r -RgC 



LxM 



f-R, 



GC 



Jh k 



RGC ~ Rn 



(6) 



In the limit 



f-R, 



GC 



large, one can write 



lim h\l ] (kp) = t- Ll - l e lkp /(kp), 



-h-likp 



p— >oo 



(7) 



where p — f— Rgc- Therefore, one obtains 



p^oo 



LV 



Y l 2 (^Roo-rJ 3h ( k \Rgc - Rn\) r^e^/ikp) . 



where 



rp(n) 

1 L'L 



(9) 
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is the T-matrix. From this, one can calculate the scattered current at distance p divided by 
the incident electron flux. In doing this, the only part of the exponential term e lkp of Eqs [7] 
and [S] to survive will be due to the imaginary part of the wave number One can write 



a e (k) = a pe (k)e~ 2 ^ 



(10) 



with the following definition for the pseudo-elastic cross section 



pe 



(*) = £lW7l*| 2 , 



(11) 



(_!)->; 



->h,h,l' 

7711,7712, — m' 



(12) 




) 



The exponential decay in Eq. (TTU]) is the result of the loss of coherence of the electron as 
it travels and is due to scatterers external to the decamer. Whereas lim^oo a e goes to zero 
for all finite £, a pe remains finite. For all practical purposes, this last quantity tends to the 
normal definition of cross section a = lim^oo a e whenever £ is much larger than the size of 
the macromolecule. 

C. R- matrix calculations 

The static-exchange approximation 4 - reduces the problem of electron scattering from a 
polyatomic molecule to a one electron problem. This approximation amounts to including 
only the ground state of the target in the close coupling expansion of the wave function, 
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and it is roughly the equivalent of the Hartree-Fock approximation 



16 



32 



45 



'or continuum states.— 



46l ; here we just sketch 



A detailed description of our method can be found in Refs. 
the main points of the treatment. 

We use the R-matrix method to solve the one electron problem.— This method consists 
in partitioning space into a short range zone, where all the channels are coupled and the 
scattering problem can in principle be treated in all its many body complexity, and an 
outer zone (external to the target electron density) in which the escaping electron only sees 
the effect of the target molecule as a multipole expansion of the electrostatic potential. In 
passing, we note that this is conceptually very suitable to the scheme we are applying in 
this work, patching together short-range scattering data for various DNA subunits, and it 
also allows us to get rid of the long-range (dipole) part of the electron molecule interactions, 
which is not relevant here. In its eigenchannel form, the R-matrix method can be formulated 
as a variational principle^ for the normal logarithmic derivative (-6) of the wavefunction on 
the reaction zone surface: 



aog(rg) _ f v W(E-H-L)*dV 

dr J v ^*S(r -r )^dV 1 ' 



where L is the Bloch operator, needed to make the Hamiltonian H Hermitian and r Q is 
the boundary between the internal and external regions. It is possible, after expanding 
the internal region wavefunction in a suitable basis set, to recast the solution of Eq. [TBI as 
an eigenvalue problem and then through basis set partitioning to shift the computational 
burden to the solution of a large linear system.— ^ As a basis set we use finite elements^ 
in all three spherical coordinates, in this way we have large but sparse matrices that are 
amenable for solution with fast sparse solvers.— 1 ^ 
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A further simplification consists in using the local density approximation (LDA) for the 
exchange potential: 

V ex (r) = ~k F F(k F ,E), (14) 
where hp is the local Fermi momentum: 

k F (r) = [37r 2 p(r)] 1/3 (15) 

and F is a functional of the energy and the local density p(r) (through the local Fermi 
momentum). The functional form we use for F is called the Hara exchange.— It has been 
extensively employed in continuum state calculations, and it is energy- dependent. The 
LDA, widely used also in density functional (DFT) calculations, 51 gives qualitatively correct 
results,— 1 ^ while it is simple enough to allow calculations for complex molecular targets. 

A polarization-correlation potential is added to this. The long range part of this potential 
is a simple multipole expansion, of which we retain only the induced dipole polarization term: 



Vpot = (16) 



where «o is the totally symmetric component of the polarizability tensor (higher order 
and anisotropic terms are much less important 46 ), and it can be calculated ah initio using 
electronic structure codes. This potential is in principle nonlocal inside the target molecule. 
We approximate it as a local potential using a form based on DFT (specifically on the 



LYP potential of Ref. l52l ) which has yielded reliable results in the work of Gianturco and 
coworkers.— This form makes use of the electron density, its gradient and Laplacian, which 
have to be calculated for each target molecule. The short and long range potentials are 
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matched unambiguously (continuously but with discontinuous derivatives) at the innermost 
crossing point, whose radius is dependent on the angles. The matching is unambiguous in 
the sense that there are two crossing points between the inner and outer potential for each 
angle, and we always choose the innermost, since the other is far in the region where the 
electron density of the molecule is very small. Choosing the outermost crossing has proven 
to give unphysical results^ in many cases. 

All the target quantities are calculated at the Hartree-Fock level using a 6-31G** basis 
set, and the target equilibrium geometries have been optimized at the same level of theory. 
The details of the calculations have been described in Ref. 16 



, including the convergence 



criteria. We notice here that the cross sections for the purines had to be recalculated due 



to an error in the calculation of the electron density, as noted in Refs. 



46 



55 



, and their 



16. These calculations are 



resonances are now shifted to lower energy with respect to Ref. 
very cumbersome, and for the level of accuracy we are aiming for here this convergence 
criterion seems adequate. 



III. RESULTS AND DISCUSSION 

A. B-form poly(A) • poly(T) base pairs decamer 

We found it judicious to start our study with the very regular poly(A) • poly(T) decamer 
(with one strand containing 10 adenines and the other 10 thymines, not considering the 
backbone) in order to avoid sequence disorder and focus more on the regularity of the spiral 
structure. This should favor comparisons with the previous toy model simulations which 
predicted important effects of internal diffraction and enhancement of capture amplitudes at 
low energy^ 1 ^. All calculations in this section were done for a coherence length of £ = 1000 
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a.u. which is roughly 15 times the decamer length. This is long enough for good coherence 
throughout the decamer and the use of a pe as a measure of the total elastic cross section. 

We used the atomic coordinates of an AT pair in the idealized structure data file ob- 
tained from the Protein Explorer internet site^ and generated the decamer by rotating this 
structure by 36 degrees and translating it by 3.38 A along the spiral axis^. Our S matrices 
were calculated using a dipole cutoff at 12 a.u., the size of our R- matrix box. The reasoning 
behind this is that the superposition of the long-range part coming from all bases would 
nearly cancel out and, in any case, be absorbed into the optical potential. The S matrices 
of adenine and thymine had to be rotated so that their principal molecular axes coincided 
with those of the corresponding bases in the decamer. Three orthogonal axes were defined 
for each base from the positions of the C2, C4, and C6 atoms of the ring part^, am for the 
R-matrix molecule positioning and for the decamer base orientation = 1,2,3). The 
following transformation is inspired by Messiah's^ treatment of rotations. We define three 
Euler angles a, /3, and 7 from the following association of 

R(a, (3, 7)^ = d Ri • a dj (17) 

with Eq. (C.45) of Messiah. The transformation matrix is then defined as 

W lmM , = e- iam r^ m ,(P)e-^ m ', (18) 

where T^ m ,(f3) is the y axis rotation matrix defined in Eq. (C.72) of Messiah. The trans- 
formed S matrix Si m y m i becomes 

§ = WSW-\ (19) 
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Let us now naively use these rotated S matrices in Eq. (J2J) using the full extent of angular 
momenta / < 8 obtained in the R-matrix calculations for an incident plane wave having a 
wave vector k perpendicular to the mid-section base-pair direction and the spiral axis. Fig. 
[TJ shows the total elastic cross section as a function of energy. One immediately notices 
the surprisingly large and suspicious values at the lower energies. The two dominant peaks 
correspond to some fifteen times the geometrical cross section of the decamer. This is not a 
low-energy dipolar effect as all dipoles are cutoff at 12 a.u. in our calculations. As we will 
see later on, the situation is even more extreme for the other decamer. The reason for this 
low-energy unruly behavior is made clear when one tries to use R-matrix results for the H 2 
molecule in solid ice^. The correct band structure of ice can only be obtained at low energy 
by cutting off the angular momenta to I < 2 instead of using the full range I < 4. This is not 
fortuitous and it can be understood by the following semi-classical argument. For an electron 
with angular momentum L ~ kr, one can write L 2 « 1(1 + 1) = k 2 r 2 . If r > d m , where 
d m is the inter-molecular distance, then obviously the electron is outside of the interaction 
space of the two molecules. This means that the relevant angular momenta are those for 
which r < d m i.e. 1(1 + l)/k 2 = 1(1 + l)/(2E e ) < d 2 m where energies are in Hartree. But 
what is the mathematical or rather the numerical reason? Let us look closer at Eq. (j2J). 
It turns out that the spherical Hankel function of the first kind h\ (kR nn >), where R nn i is 
the distance between the two scatterers n and n', diverges as (/c.R nn ')~^ 1+1 ) for small values 
of the argument (kR nn /). The singularity is dominated by the nearest neighbor distance 
d m . For /, V 2 = 4, li can be as large as 21 = 8. The reason for this comes from the angular 
momentum composition rule that says the sum over l\ can reach the uppermost value of 
l\ < \l + l' 2 \. So even though one expects the V 2 = 4 components of the T- matrix T = S — 1 
to get smaller as the energy decreases, the product of h\^(kR nn >) with T gets overly large 
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because of the much more singular behavior of the Hankel function. The situation gets much 
more acute in our case since for I < 8 the sum over l\ goes up to l\ = 16. 

One thus also needs to use a cutoff to describe multiple collisions in the decamer. Only 
values I < l Q should be retained such that E(l a , d m ) < E e < E(l Q + 1, d m ) where E(l Q , d m ) = 
1(1 + l)/(2d^ n ). One might think that d m should be equal to the base stacking distance along 
the spiral axis direction. But it is clear that an electron that scatters from one molecule to 
another does so from every part of the first molecular subunit to every part of the other, 
including from one end of the one molecule to the diagonally opposite end of the other one. 
We have thus chosen, with some arbitrariness, the value d m — 11 a.u. which is of the order 
of this distance, the size of the bases, the distance between base centers in the base pairs and 
close enough to the size of the R-matrix sphere to retain all of its important energy dependent 
characteristics. But this now poses a new problem having to do with the discreteness of /. A 
strict imposition of these cutoffs would result in a piecewise chopped cross section as can be 
seen in Fig. [TJ One way to circumvent this difficulty, which we have adopted throughout, is 
to interpolate any scalar quantity, calculated for all integer values of angular momenta, at the 
non- integer value of I obtained from the solution of E e = 1(1 + 1)/ (2d^ n ) . A similar procedure 
is being used in a study of a H 2 dimer— and the resulting elastic cross section is found to 
be within 5% of the one calculated using the R-matrix for energies larger than 2.5 eV. The 
result of this procedure on the decamer is shown in Fig. [2j A comparison with the bare 
molecular cross sections obtained from the R-matrix calculation shows good concordance of 
most of the peaks, with the exception of the one at 0.24 H in thymine which is lessened and 
perhaps split. This is somewhat reassuring for the credibility of the interpolation procedure. 



Let us now look for signs of internal diffraction. As discussed in Refs 



28 and 



29, this 



would occur because of the regular spacing of base pairs along the spiral axis. We have 
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calculated various quantities which were previously shown to be sensitive to this regular 
spacing. Fig. [3](a) shows the quantity T w (0), equal to the average of the square of the 
electron wave function on the bases. It is compared to the one obtained by dilating the 
distances along the spiral axis by a factor 1.05. One expects a down shift in energy for 
those peaks sensitive to the inter-base-pair spacing d\\ by some 10 % (E oc k 2 oc A~ 2 oc d^ 2 ). 
Three peaks seem to behave this way although the one at 0.42 H at regular spacing is more 
probably the internal diffraction one that was deduced in the previous publications. The 
unshifted peak slightly below 0.1 H is of local origin. Fig. H]shows a contour plot of 7(0, R n ), 
the absolute square of the wave function at each base's center, in energy and base sequence 
space. One must of course realize that these contour lines are calculated by interpolating 
7(0, R n ), which is defined only for a discrete set of the index n. The virtue of such a graph 
is to enable a quick survey of cool or hot spots at which the square of the wave function 
is small or large. The sequencing index varies from 1 to 10 going up the A strand and 
from 11 to 20 going down the T strand. The peak structures in r u) (0) showing appreciable 
enhancement of the wave function are seen to occur mostly on the base pairs that lie nearly 
perpendicular to the incident electron direction (actually, the closest to perpendicular is at 
an angle of 90 ± 18 degrees), at the beginning, middle, and end of the decamer, where phase 
coherence is more favorable. We have also calculated the axially scattered current per unit 
solid angle in the +z direction Fig. [5]^a) and the forward scattered current per unit solid 
angle in Fig. 0(b) under identical axial spacing conditions as in Fig. [3](a). The only shifting 
peak that is common to all figures is the one at 0.42 H at regular spacing. This is clearly 
the internal diffraction peak. One should note the large peak at 0.34 H in Fig. 0(a). There 
is nothing special that can be seen in the wave function in Fig. H] at that energy. This result 
can only be understood by interference between the scattered beams emanating from the 
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bases. For the current along the axis, aside from the contribution of the beam amplitudes 
at each base R n , there is an extra exit phase factor at position z — > oo on the axis outside 
the decamer, which is proportional to exp[ik(z — z n )]. This means that there is an optimal 
value of k at which exp[ik(z n ' — z n )] between neighboring base pairs along the spiral axis 
will "synchronize" the beams and produce an overall good phase coherence. The position 
of the peak yields k ~ 0.825 and a A ~ 7.6. The vertical distance between base pairs is 
6.4 a.u.. Why don't the numbers match? This is because we are dealing with a complex 
superposition of many different partial waves with different phases for all I (and rn) less or 
equal to 8. The optimal value of A has to compromise with all of these and this happens for 
A ~ 7.6. 

We can now also answer the question whether or not there is appreciable enhancement 
of the partial capture factors at low energy. Fig. [3(b) clearly says so for energies less than 
say 0.1 H, a region of interest for the low-energy I = 2, 3 shape resonances of the bases (see 
Fig Eta) showing cross sections of the 4 bases). 

B. B-form GCGAATTGGC base pairs decamer 

Now that we know how to look for internal diffraction, let us repeat the analysis of the 
last section on our sequence disordered decamer for £ = 1000 a.u. and k perpendicular to 
the mid-section base-pair direction and the spiral axis. The total elastic cross section as a 
function of energy using the full extent of angular momenta I < 8 exhibits a cross section 
that has peaks some 500 times the geometrical one. This confirms the unphysical behavior of 
this procedure at low energy. Using the interpolation procedure, however, restores a credible 
cross section as shown in Fig. E(a). A comparison with the cross section of individual bases 
again reveals a satisfactory correspondence of most peak structures. Fig. [7] shows the 
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quantity T u ,(0), equal to the average of the square of the electron wave function on the 
bases. It is compared to the one obtained by dilating the distances along the spiral axis by 
a factor 1.05. The diffraction signature around 0.4 H is still present but with an appreciably 
reduced enhancement factor relative to the plane wave value of 1. This is the effect of 
sequence disorder which turns out to be stronger than what was anticipated from the toy 
model calculations. The lower energy structures are similar to those of the poly(A) ■ poly(T) 
decamer. Fig. E](a) shows a contour plot of 7(0, R n ), the absolute square of the wave function 
at each base's center, in energy and base sequence space. The sequencing index varies from 1 
to 10 going up the GCGAATTGGC strand and from 11 to 20 going down the complementary 
strand. The peak structures in r„,(0) of Fig. [7] showing appreciable enhancement of the 
wave function are seen to occur mostly in the mid region, on the base pairs that lie mostly 
perpendicular to the incident electron direction. Sequence disorder has a visible effect on the 
interference patterns within the decamer and not only on internal diffraction. It somewhat 
desynchronizes the ends from the mid-section. The somewhat irregular distribution of the 
peaks at lower energies (E < 0.3 H) in both decamers studied thus far shows that these 
structures are local ones and perhaps a signature of weak localization. Localization has 
been predicted in poly(dA)-poly(dT) decamer studies due to structural changes promoted by 
thermal fluctuations 365 . There is a parallel to be made between the time evolution of localized 
states and our energy dependence. The axially scattered current in the +z direction (from 
the G end to the C end) shown in Fig. [H](a) also presents a clear but quite diminished internal 
diffraction peak. But the large peak at 0.34 H in Fig. 0(a) has now disappeared, stressing 
the importance of sequence disorder. The forward scattered current of Fig. Mjo) shows no 
evidence of internal diffraction. This parameter is less sensitive to internal diffraction as was 
already evident from the poly(A) ■ poly(T) results. The weighted partial capture factors are 
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very similar to the the poly(A) • poly(T) one. 

Let us now reduce the coherence length to £ = 20 a.u. to simulate phase decoherence 
within the decamer. The pseudo cross section a pe of Eq. (TTOT) is shown in Fig. Mjo). It 
is very jagged compared to the previous cross sections. Remember that the values are 
deceptive since a e decreases by a factor of fifty at a distance of 40 a.u. from the geometrical 
center of the decamer and only a few a.u. from the ends. A similar behavior for T W (L) can 
be understood in terms of huge localization resonances very apparent in Fig. ffljo), which 
shows the contour plot of 7(0, R n ), the absolute square of the wave function at each base's 
center, in energy and base sequence space. Such localization resonances were seen in highly 
disordered single- strands of DNA— although not as pronounced. There is localization at 
all energies although it gets weaker as the energy increases. There is no obvious pattern in 
their position within the decamer. The enhancement of the partial capture factors is very 
large in all resonances with energy less than 0.3 H. These large numbers are an artifact 
of a large imaginary part to the wave number which prevents bases from "seeing" their 
far neighbors and thus forces a network of strong correlation of amplitudes with the near 
neighbors. This is understandable in the presence of important inelastic losses. This appears 
inappropriate, however, for disorder due to stray elastic scatterers since correlations will still 
exist with all bases even though there is randomization of the local phases. These correlations 
should weaken this network. It would be important to verify this by adding some genuine 
randomness such as including water molecules bound to DNA. We plan to do this in the 
near future. 
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C. A-form GCGAATTGGC base pairs decamer 

We finally look at the A-form of the GCGAATTGGC base pairs decamer. There is 
considerable tilt (22.6 degree inclination) in addition to roll (-10.5 degrees propeller twist) 
of the base pairs in this formal Moreover, the pairs are more tightly packed having a rise of 
2.54 A compared to 3.38 A for the B-form and there are 11 base pairs per turn of the helix. 
It is of interest to look for the differences in the scattering results. 

Fig. [TD] shows the cross section for this A-form decamer. One immediately notices the 
globally smaller values compared to its B-form analog in Fig. 0(a). This was unexpected in 
view of the closer packing of the base pairs and of the results from the toy modeP^ which led 
us to expect larger capture factors. This emphasizes the need for realistic scattering models. 
It would seem that there is some increasing destructive interference as the base pairs get 
closer. In order to validate this hypothesis, we have unwound the A-form by increasing the 
rise and helical twist to agree with the B-form values, but keeping the base pair tilt and roll 
intact. Fig. [10] shows the cross section for this unwound decamer. One recovers the same 
scale of values as for the B-form decamer. It is therefore indeed destructive interference 
caused by the closer proximity of the base pairs that is responsible for the smaller cross 
section of the A-form and not the tilt and roll. The latter are, however, surely responsible 
for the difference in peak structures. Note that there is no internal diffraction to be seen in 
the A-form in the chosen energy range. The distance between nearest base pairs (the rise) 
is too small which pushes the diffraction effects to energies beyond 0.5 H. 

Finally, we have computed scattering parameters for both the A and B forms which are 
more easily comparable with existing experimental data. The averaged electron current 
scattered along the DNA axis (±Z scattered directions and X, Y incident directions) is 
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shown in Fig. [TTJ as a function of electron energy. Similarly, the magnitude of the square 
of the wavefunction averaged over all bases is shown in Fig. [12j To facilitate comparison 
with experimental data, the energy scale has been converted to eV and the zero electron 
energy shifted by the estimated value U op of the polarization energy. This shift is necessary, 
since in thin film experiments with electrons incoming from vacuum, the zero energy is 
conveniently referenced to that of the vacuum level. The U op values for the A and B forms 
were evaluated by estimating the polarization energy on a base due to its near neighbors. 
As we wish to make comparisons with experimental results on DNA, we have added an 
interaction with the nearest phosphate and deoxyribose subunits as well as to some four- 
close structural water molecules that would be present in real DNA. We have utilized the 
formula U op ~ — ^a/(2d 4 ) where a was given an average value of 100 (a.u.) 3 for the in- 
plane polarizability of the partner base^ -, an average value of 47 (a.u.) 3 for the out of plane 
polarizability of the two axial neighbors- , an average value of 5 and 80 (a.u.) 3 for the the 
isotropic polarizabilities of phosphate and deoxyribose respectively^, and of 10 (a.u.) 3 for 
the water molecules^. We used a distance d ~ 11 a.u. for the partner distance and 7 a.u. 
for all other distances^ except for the A-form axial nearest neighbor distance, which was 
taken to be 5.3 a.u.. Needless to say that the values obtained for U op , -1.5 eV and -2.5 eV 
for the B-form and A-form respectively, are approximate although credible. 

IV. COMPARISON BETWEEN THEORY AND EXPERIMENTS 

Low-energy electron (LEE) experiments with DNA have been performed on dry films, 
in ultra-high vacuum, where the molecule adopts the A-form owing to non-structural water 
loss^. These experiments have measured the damage inflicted to DNA mostly by electrons 
in the 0-20 eV range, i.e., the yield functions, in the form of base^ 1 ^, sugar— and phosphate 
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lesions^ 6 -, single and double strand breaks^ 7 - 1 ^ and base release^. Thus, no electron scat- 
tering experiments are presently available to be directly compared to our theoretical results. 
However, some of the mechanisms, which have been invoked to account for the magnitude 
of the yield of specific damages and peak energies in their yield functions below 15 eV, rely 
on scattering properties of LEE within the DNA molecule. 

In experiments with thin films of plasmid DNA the yield function for single strand breaks 
(ssb) exhibits maxima at 0.8, 2.2 and 10 eV, with a shoulder at 6 eV—^^ which appears 
as a distinctive peak at this energy in the ssb yield function of synthetic single stranded 
DNA films 70 . The peak positions were determined with an accuracy of 0.3 eV. From the 
analysis of LEE-induced products from this latter type of films, the yield function for base 
release was also found to exhibit peaks at 6 and 10 eV— . Below 15 eV, the yield function 
for the induction of double strand breaks (dsb) was found to be dominated by a single 
peak located at 10 eV- 1 ^. Owing to their low energies the 0.8 and 2.2 eV maxima could 
easily be interpreted as shape resonances^, but the exact mechanism leading to rupture of 
the C-0 phosphodiester bond (i.e., ssb) 6 - is not obvious for two reasons. First, the cross 
section for ssb induction below 3 eV is of the order of 10-17 cm2 per base pair—, which 
is fairly large for damage caused by a single anion dissociative state at 0.8 and 2.2 eV, 
respectively. In fact, this value is almost as large as that measured at 100 eV, in the same 
DNA, where a plethora of ionization and dissociation channels are available 6 ^ 1 ^. Coherent 
enhancement of the electron wavefunction within DNA may therefore increase at low energies 
the electron capture cross section of transient anion states responsible for ssb. Secondly, in 
the experiment of Martin et al.— , the 0.8 and 2.2 eV peaks did not coincide with the energy 
of the dissociative phosphate anion which is known to rupture the C-0 bond of the DNA 
backbone 6 ^.. Instead, these peaks coincided with those in the electron capture cross section 

23 



of the DNA bases. The latter result lead Martin et al.— to postulate that the electron was 
first captured by the basis and then transferred to the phosphate group. 

In fact, according to the theoretical studies of Simons' group, below 3 eV electrons 
cleave the CO bond of the DNA backbone at the 3' and 5' positions essentially via electron 
transfer-^i£^i. An electron is first captured by a base to form a pi* transient anion and 
afterwards the additional electron transfers, via the sugar moiety, to the phosphate unit 
where it occupies a sigma* orbital at either the 3' or 5' C-0 positions. The resulting anion 
state being dissociative leads to C-0 bond cleavage (i.e., a ssb). Electron transfer is not 
necessarily limited to the nucleotide where capture occurs^ 5 - 1 ^. Thus, if the main mechanism 
leading to ssb is electron transfer from the bases, the large cross section for damage below 
3 eV could be explained by invoking strong constructive interference of the electron wave 
function due to stacking of the bases along the DNA chain. 

This cross section could maximize at 0.8 and 2.2 eV, if the incoming electron preferably 
scatters along the DNA axis at these energies, where it would cause preferential coherent 
enhancement due to base stacking. This is exactly what is found in the energy dependence 
of the magnitude of the electron current scattered along the axis of A-DNA in Fig. [TTJ 
the current in the Z-direction is maximized at 0.8 and 2.2 eV. When DNA is modified to 
its B form the two peaks move closer to each other and little correlation is found with 
the experimental values. In Fig. [121 the average of the square of the scattered electron 
wavefunction on a base is found to maximize at 1.3 and 2.8 eV in fair agreement with the 
experimental maxima in the yield function for ssb. For the B form, Fig. [TJ] shows no 
two-peak correlation between theory and experiment. Thus, by including both the shape 
resonance wavefunctions and constructive interference due to base stacking our calculation 
can represent fairly well the energy dependence of the yield of ssb below 3 eV. 
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The broad 6-eV feature in the experimental ssb yield function of DNA, which spreads 
from 5 to 7 eV— , is also in good agreement with the theoretical results of Figs [TT1 and [121 for 
the A-form. Both curves exhibit two peaks around 6 eV, which if unresolved would produce 
a broad maximum around 6 eV, as observed experimentally. 

The 6-eV feature has been studied in detail by Zheng et al. who bombarded thin molec- 
ular films of a short single strand of DNA, with electrons of energies between 4 and 15 
eV— . By high-pressure liquid chromatography, they identified 12 fragments of the oligonu- 
cleotide GCAT sequentially composed of the bases guanine (G), cytosine (C), adenine (A) 
and thymine (T). The yield functions exhibited maxima at 6 and 10-12 eV, which were 
interpreted as due to the formation of transient anions leading to fragmentation. Later, 
they analyzed the products induced by 4-15 eV electrons incident on two abasic forms of 
the tetramer GCAT, i.e., XCAT and GCXT, where X represents the base, which has been 
removed and replaced by a hydrogen atom . The results obtained at an incident energy of 6 
eV showed that essentially no strand break occurs at positions in the backbone corresponding 
to those of the missing base. This finding clearly indicated that at 6 eV, and possibly below, 
electrons break the DNA backbone almost exclusively via electron transfer. Furthermore, the 
total yield of all the bases released and ssb induced by electrons were found to be strongly 
affected by the presence of an abasic site; in both XCAT and GCXT, the yield of detached 
bases was found to be up to an order of magnitude smaller than that from GCAT 77 . Thus, 
the initial electron capture amplitude was suggested to be highly sensitive to the number 
and possibly the geometrical arrangement of the bases, indicating the presence of a strong 
collective effect. According to the results of Figs [11] and [121 this collective effect could be 
related to strong electron scattering on axis around 6 eV, which exhibits an 8-fold increase 
in magnitude in going from 4 to 6 eV. Coherence enhancement of the scattered electron 
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wavefunction around 6 eV could also play an important role. Both phenomena should be 
highly sensitive to base removal, since they are directly related to base stacking and their 
periodicity. 

Finally, the maximum at 10 eV in the electron energy dependence of the calculated 
averaged axial scattered current and square of the wavefunction averaged on all bases for the 
A-form correlate well with the strong maximum found experimentally in the yield function 
for ssb and dsb induced by LEE impact on plasmid DNA- 1 ^. This result suggests that 
coherence also influences DNA damage in the 10 eV region; i.e., at energies where the 
formation of core excited resonances is the dominant mechanism implicated in bond rupture^. 
Such anion states, which consist of two electrons occupying electronically excited orbitals 
around a positive hole of a DNA subunit, are highly localized^. It is not obvious how such 
anion states would couple to a diffracted wave along the DNA axis, a problem which has 
not been addressed in the present work. Furthermore, present and previous^ 1 ^ partial wave 
analysis of the enhanced electron capture probability on specific subunits due to diffraction 
has been show to be much reduced in the 10 eV region compared to low energies. It is 
therefore possible that the energy coincidence mentioned above is fortuitous. 

To check the validity of all of our comparisons, we calculated the averaged electron 
current scattered along the DNA axis (Z direction) and the magnitude of the square of the 
wavefunction averaged over all bases for different base sequences of the A-form. We observe 
a fairly good stability in the features of the curves up to about 7 eV. In other words, the 
effect of diffraction is not very dependent on sequence and therefore a comparison with 
experimental data obtained with different sequences can be considered significant below 
7 eV. Thus, the stability of the structural information in Figs [11] and [12] with different 
base sequences and the comparison with experimental results indicates that wave function 
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interference should be taken into account to describe the mechanism of action of electrons 
with energies lower than 7 eV in DNA. Beyond this value, the energy of the calculated 
minima and maxima change according to sequence. Thus in this case, multiple scattering of 
the electron wave in DNA is highly sequence dependent and comparison with data obtained 
with plasmid DNA, which is longer and has a different sequence, is not considered to be 
significant. 



V. CONCLUSIONS 

This article describes the first attempt to use high-quality gas phase scattering data to 
gain insight into the interaction, in condensed phase, of DNA and low energy electrons, 
which is relevant for radiation damage to nucleic acids by ionizing radiation. In particular, 
we used R-matrix calculations performed by some of us^ for the gas phase scattering. To 
explore the interaction of a continuum electron with the DNA double strand, we couple the 
gas phase calculations with a multiple scattering framework developed by some of ua 2§i2iL22j2l 
a few years ago. The combination of R-matrix calculations and multiple scattering has been 
recently shown^^ to give accurate results for simpler systems as the water dimer and ice. 
These recent works have also inspired the angular momentum cutoff procedure and the 
interpolation of the scalar quantities described in Sec. IIHI 
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show that some peaks 



The results we show here, in a similar fashion to Refs. 
in electron scattering from macromolecules can be due to diffraction instead of formation 
of negative anions. They also show that single basis calculations are not enough to gain 
an understanding of the process and that collective effects can be as important as local 
ones. Predictably, disorder is shown to play an important role, as evidenced from the effect 
of sequence disorder and of electron decoherence that we have performed introducing an 
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imaginary part of the optical potential, as in Sec. flTA] It is difficult to give a quantitative 
assessment of the role of disorder in electron scattering from biomolecules within our current 
model, but it would be definitely an interesting topic unto itself, also for investigations 
by the experimental community, maybe performing measurements in solution at different 
temperatures. We have shown how decoherence leads to strong resonance localization and 
large wave function amplitudes, generating strong resonance enhancements across the energy 
range we explored. 

Despite the simplicity of our model, we found fairly good correlation between calculated 
values for the A form and experimental results obtained on LEE-induced damage to DNA 
in this configuration. Below 7 eV, the calculated averaged electron current scattered along 
the DNA axis (Z direction) and the magnitude of the square of the electron wavefunction, 
averaged over all bases for different base sequences of the A-form, were found to be essentially 
independent of the nature and sequence of the bases. It was therefore possible to compare 
calculated values with the results of experiments performed with a longer DNA molecule of a 
different sequence. The structures which appeared in the calculated current scattered along 
the DNA axis and the square of the wavefunction correlated well with the maxima found 
in the yield function for breaking a single and two adjacent strands of DNA by the impact 
of 0-7 eV electrons. The correlation indicates that a substantial increase in DNA damage 
occurs at the energies of preferential scattering of the incoming electron along the DNA 
axis; i.e., at energies where increased coherent enhancement occurs due to base stacking. 
Although our results for A-DNA and B-DNA are not markedly different, theory-experiment 
correlation appears to be better for the A form, which is likely to be the one present in thin 
film experiments. 29 It appears from both comparisons, that diffraction effects within DNA 
should be taken into account to describe the mechanism of action of electrons with energies 
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below 7 eV. This may also hold true for higher energies, but it is difficult to compare the 
present calculations with experimental data beyond 7 eV, because core-excited resonances 
implicated in the damaging process were not included in our model and because the energy 
dependence of the scattered current and of the square of the wavefunction change according 
to sequence. 

Our model constitutes an attempt to relate the gas phase and condensed phase areas of 
the research in electron scattering from biomolecules. The model is far from perfect since 
it has many assumptions. In particular it is difficult to draw solid conclusions on radiation 
damage from this model, first of all because we consider only elastic scattering events, and 
completely neglect the motion of the nuclei, which would add a prohibitive new level of 
complexity. Also, considering DNA as rigid is far too simplistic in solution but also to 
some extent in thin films where at best the material can be considered "amorphous" and 
structural water is present anyway. The motion of the electron in the liquid should also be 
considered, and it would be best to use an electron distribution taken from photoemission in 
a biological medium like water instead of using a plane wave^. One simple improvement to 
the model could be to consider structural water, with positions taken possibly from relevant 
DNA crystal structures, and we are planning to explore this possibility to gain a better 
understanding of the role of disorder and parasite scatterers. 
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FIG. 1: Total elastic cross section of the poly(A) -poly(T) decamer as a function of incident electron 
energy using the full angular momentum content of the scattering matrix (full line) or imposing a 
strict application of the angular momentum cutoff criterion (dark circles). 
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FIG. 2: Interpolated total elastic cross section of the poly(A) -poly(T) decamer compared with the 
single adenine and thymine R-matrix cross section values as a function of incident electron energy. 
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FIG. 3: Interpolated base averaged value of the square of the electron wave function for the normal 
spacing and with an expanded spacing 1.05 times the regular one (a) and interpolated value of the 
base averaged weighted partial capture factors (b) for the poly(A) • poly(T) decamer as a function 
of incident electron energy. 
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FIG. 4: Contour plot of the square of the wave function at each base's center for the poly(A)-poly(T) 
decamer as a function of energy and base sequence number n. Contour lines are equally spaced at 
integer values. 
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FIG. 5: Interpolated value of the axial scattered current in the +z direction (a) and of the forward 
scattered current (b) of the poly(A) • poly(T) decamer as a function of incident electron energy for 
the normal spacing and with an expanded spacing 1.05 times the regular one. 




0.2 0.4 

Energy (Hartree) 



0.2 0.4 

Energy (Hartree) 



FIG. 6: Interpolated total elastic cross section of the GCGAATTGGC B-form decamer at £ = 1000 
compared with the single adenine, cytosine, guanine, and thymine R-matrix cross section values 
(a) and at £ = 20 (b) as a function of incident electron energy. 
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FIG. 7: Interpolated value of the average of the square of the electron wave function on the bases 
of the GCGAATTGGC B-form decamer at £ = 1000 as a function of incident electron energy for 
the normal spacing and with an expanded spacing 1.05 times the regular one. 
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FIG. 8: Contour plot of the square of the wave function at each base's center for the GC- 
GAATTGGC B-form decamer at £ = 1000 (a) and at £ = 20 (b) as a function of energy and 
base sequence number n. Contour lines are equally spaced at integer values in (a) and multiples 
of 5 in (b) 
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FIG. 9: Interpolated value of the axial scattered current in the +z direction (a) and the forward 
scattered current (b) of the GCGAATTGGC B-form decamer at £ = 1000 as a function of incident 
electron energy for the normal spacing and with an expanded spacing 1.05 times the regular one. 
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FIG. 10: Interpolated total elastic cross section of the GCGAATTGGC A-form decamer at £ = 
1000 (open circles) and of the unwound counterpart having B-form rise and helical twist (black 
circles) as a function of incident electron energy. 
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FIG. 11: Average of the scattered current in the two axial directions and for two incident and 
mutually orthogonal directions of the GCGAATTGGC A-form and B-form decamers at £ = 1000 
as a function of incident electron energy. 
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FIG. 12: Average of the square of the electron wave function on the bases over two incident and 
mutually orthogonal directions of the GCGAATTGGC A-form and B-form decamers at £ = 1000 
as a function of incident electron energy. 
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